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Amplification induced by a periodic coefficient and a multiplicative white noise term 
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The author studied the amplification of the the amplitude in a Mathieu-like equation with multi- 
plicative white noise. An expression of the exponent on parametrically amplified regions was derived 
theoretically by introducing the width of time interval, and the exponents were calculated numeri- 
cally by solving the stochastic differential equations by a symplectic numerical method. It was shown 
that the exponent decreases with a parameter a, reaches the minimum and increases after that, where 
the value of a is determined by the intensity of noise and the strength of the coupling between the 
variable and the noise. The author found that the exponent as a function of a has only one minimum 
at a + on parametrically amplified regions of a — 0. This minimum value is obtained theoret- 
ically and numerically. The existence of the minimum at a + indicates the suppression of the 
amplification by white noise. 
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1. Introduction 

In past few decades, many researchers have investigated the roles of noise, and then remarkable 
phenomena were found. For examples, such phenomena are stochastic resonance, ! ~ 4) phase transition 
induced by multiplicative noise 5) and so on. 6_10) A basic system in which multiplicative noise acts 
is an oscillator with varying mass and/or friction. The amplification by noise is important, because it 
brings the unstableness of the state, the clear signal and so on. 

The role of noise acting on a harmonic oscillator multiplicatively was investigated by 
Stratonovich. n) He found that the amplitude of the oscillator with friction is amplified by noise. 
Recently, the phase transition was investigated for pendulum with a randomly vibrating suspension 
axis. 12 ' Mallick and Marcq investigated nonlinear oscillators with white or colored noise. 13 ~ 15) They 
showed that the time dependence of the amplification obeys power law. 

Another mechanism of the amplification is parametric resonance. 16) Some physicists have already 
been investigated the effects of additive white noise acting on the harmonic oscillator with a periodic 
coefficient. 17) They showed that the covariance of a coordinate has a minimum value at a nonzero value 
of the parameter which is included in the periodic coefficient. In the study of a parametric oscillator 
driven by a periodic coefficient with additive white noise, 18 ' 19) the suppression of the fluctuation was 
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found. The parametric resonance induced by multiplicative colored noise was investigated in ref. 20. 
In field theories, it was shown that a field is amplified by another oscillating field. The presence of 
noise leads to an increase in the rate of particle production. 21 ^ In <p A theory, the equation of motion of 
the field of a nonzero mode near the bottom of the effective potential may have an oscillating mass 
term with noise. 22) In such a system, the magnitude of the amplification is modified by white noise 
quantitatively. A Mathieu-like equation which has a multiplicative noise term may appear even in 
the end of phase transition, as shown in ref. 22. Experimentally, some physical systems which are 
described by the equations with a periodic coefficient and a multiplicative noise term were studied. 23) 
The equations to describe such systems are similar to the equation of pendulum. 

The differential equation with a periodic coefficient and a multiplicative noise term appears in 
many physical systems. Multiplicative noise may amplify or suppress the amplitude, as found in the 
case of additive noise. The effects of multiplicative white noise on the system which shows the para- 
metric amplification should be investigated, because the suppression of the amplitude is directly related 
to the stability of the system. 

In this paper, I studied the effects of multiplicative white noise acting on the system that the para- 
metric resonance occurs when the values of the parameters of the system are adequate. A stochastic 
differential equation was analyzed by introducing the width of time interval, and I derived an approx- 
imate expression of the exponent on parametrically amplified regions where parametric resonance oc- 
curs when no noise exists. The behavior of the exponent was discussed qualitatively and the minimum 
value of the exponent as a function of the parameter a was estimated, where this parameter indicates 
the intensity of noise and the strength of the coupling between the variable and the noise. The stochas- 
tic differential equations were solved numerically by a symplectic method to avoid the amplification 
by numerical error, and the exponent was extracted from the average of the trajectories. The behavior 
of the exponent as a function of a was displayed numerically with the other parameters fixed. I found 
that the exponent has only one minimum at a + on parametrically amplified regions of a - 0. The 
existence of the minimum indicates the suppression of the amplification by white noise. The results 
provide insight in the systems with periodically varying parameters and multiplicative noise. 

2. Qualitative properties of the exponent on parametric amplified region 

2. 1 The basic equation 

The equations with a periodic coefficient and a multiplicative white noise term are interested in 
many branches of physics. We come up with the following equation which describes some phenomena 
(under some conditions), for examples, the motion of surface wave, 23) the time variation of the charge 
in a capacitor, 23) reheating in inflationary cosmology, 21) motion of the condensate near the potential 
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well 22) and so on: 

^- + [co 2 + B cos (Qt + 9) + Dn(t)] = 0, (1) 

where to, B, D, Q and 8 are constants, and t is time. The quantity n(f) is white noise with the property, 
(n(t)) - and (n(t)n(t')) - 8(t—t), where the notation (• • • ) represents statistical average. By applying 
the transformation z = (ot + Goj/Q., this equation is transformed to 

d 2 d> r 

+ [ 1 + B cos (yz) + ar{z)\ = 0, (2) 
where the quantities a, B and y are defined as follows: 

a = D/aj 3/2 , B = B/co 2 , y = Q/co. (3) 
The quantity r(z) is defined by r(z) = to~ l/2 n(z/co - Q/Q). Therefore r(z) has the following properties: 

<Kz)> = 0, (r(z)r(z')) = 6(z-z'). (4) 

The starting point in this paper is eq. (2) with eq. (4). 

Equation (2) is rewritten with the variable which is defined by - dcp/dz: 

d(p = p^dz, (5a) 

dp<t, = - [1 +jScos (yz)] (pdz - a<p o JW. (5b) 

The quantity W(z) is defined by W(z) - f ds r(s) and this is a wiener process, where the quantity 
Zo is an initial time. (Here, the symbol o represents Stratonovich product.) I attempt to solve eqs. (5a) 
and (5b) numerically in § 3. 

Equation (2) is just a Mathieu equation when a is zero, and this equation has resonance bands. 
The Mathieu equation corresponding to eq. (2) with the relation 2u - yz is given by 

d 2 d> 

—r + (a-2q cos(2w)) = 0, (6) 
du 1 

where a = 4/y 2 and -2q = AB/y 2 . Then the bands are distinguished by positive integer n with the 
relation n 2 = 4/y 2 . Therefore the values of y in resonance bands are close to 2/n at a - 0. 



2.2 Minimum Value of the Exponent on parametrically amplified regions 

In this subsection, I attempt to estimate the growth rate of the field amplitude in time. I call this 
rate the exponent in this paper. This rate is given by limsup z~ X In [|(0(z))| / |0(O)|] , where 0(0) is the 
initial value. I use the solution of the Mathieu equation in order to solve eq. (2) approximately in the 
region where parametric resonance occurs at a - 0. The equation at a = is 

<I> + [l+/?cos(yz)]<D = 0, (7) 
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where the dot represents the derivative with respect to z. The quantity (p is represented as a product of 
CD multiplied by a new variable if/: 

<f> = ®if/. (8) 

By substituting eq. (8) into eq. (2), if/ satisfies the subsequent equation: 

if/ + 2 (o/o) iff + ar(z)if/ = 0. (9) 

The exponent of CD was investigated by many researchers in detail. Therefore, the exponent of <f> is 
estimated by obtaining the exponent of ip approximately. 

Here I denote the exponent of <p at a = as s = s(J3,y) which is just the exponent of CD. The 
time dependence of O is obtained by solving eq. (7). One method to solve approximately in the first 
resonance band is performed by putting the form of CD with the assumption P ~ as follows: 16 ' 24 ~ 26) 

0) = P{z)e iyz/2 + P*{z)e- iyzl2 . (10) 

This solution gives the approximate expression of CD. From this calculation, I obtain 

G> ~ e^F^z), (11a) 

Fi(z) := Cj yzl1 + C*e~ iyzl2 , (lib) 

where C is a complex constant and s\ is the exponent. It is conjectured that the exponent s\ is close to 
the exponent s in the first resonance band. With eqs. (11a) and (lib), I obtain 

cb/cD ~ si + P\/F\. (12) 

The exponent is estimated by solving eq. (9) with eq. (12). However, it is not easy to handle eq. (9). 
Instead, in eq. (9), I replace (6/cl>) by the average of (0/0) in time. This replacement is probably 
valid when the change of if/ is slow enough. The condition that this replacement is valid is extracted as 
follows. The equation for if/ is described as 
d_ 

Tz 



t d 
(F,) 2 —ip 
dz 



+ 2s, 



i d 
(Fx) 2 —if/ 

dz 



+ ar(z)(F l ) 2 if/ = 0. (13) 



The solution for small a will be close to the solution at a = which satisfies (F\) 2 if/ = C exp (-2i , 1 z), 
where C is constant. Here I denote the period of F\(z) as T and the time scale of the exponential 
decrease as T*. The variaton of if/ is slow when the absolute value of ip is small. It is obvious from this 
approximate solution, (Fi) if/ - C exp (-2siz), that the absolute value of if/ is small when the function 
F\(z) is approximately constant in the width T*. Therefore the solution obtained by this replacement 
will be valid when the condition, T » T* , is satisfied. I put F\(z) = cos(yz/2) in order to check this 
condition. In this case, T is Any' 1 . Then this condition is expressed by the inequality 8ns » y, because 
T* is estimated as 1 /(2s). The value of y on the first resonance band is close to 2 and the value of si is 
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about 0.45 at ft = 2 and y - 2. Then this inequality holds at a = 0. Therefore the replacement by the 
average should be available to extract the properties of the exponent for small a. Keeping these facts 
in mind, I estimate the averaged value of the exponent iff. The average of O/O in one period of F\(z) is 
equal to si. Therefore, the approximate equation for iff under this approximation in the first resonance 
band is 



iff + 2s \iff + ar(z)iff = 0. 



(14) 



I derive an approximate equation in other resonance bands in the similar way. The variable O is 
expanded as 



O = ^ [P n (z)e in ^ 2 + P* n (z)e- in ? z/2 ] + R(z). 



(15) 



n=\ 



The growth of the function P m (z) is largest in mth resonance bands. Therefore, O in the mth band is 
approximately given by 

<D ~ e s >» z F m (z), F m (z) := Ce im ^ 2 + C *e~ im ^ 2 , (16) 
where s m is the exponent. The equation for iff becomes 

iff + 2s m iff + ar(z)iff = 0. (17) 
Then I attempt to solve the subsequent equation for if/ in the resonance band: 

iff + 2sip + ar(z)if/ - 0. (18) 
The purpose in this paper is to estimate the exponent. Then I put the form of if/ as follows: 

iff = iff exp^J dz'<r{z')\. (19) 



Substituting eq. (19) into eq. (18), I obtain the equation for cr. 

& + cr 2 + 2scr + ar(z) = 0. 



(20) 



At first, I find the solution when r(z) is constant. The solution of eq. (20) is categorized by the 
quantity D which is defined as As 2 - 4ar. For D > 0, 1 have 



P dz'a-(z') = (s + Vo/2) (z - zo) + 

•JZQ 



1 



VI) 



In 



1 -Ce 



1 -Ce~^< 



(21) 



where C is an arbitrary constant. For V5)z » 1, the last term of the right-hand side in eq. (21) is 
negligible for the exponent. For D < 0, 1 have 

cos(^zo + C') 



f 



dz'cr(z) = -s (z - Zo) + In 



cos 



(22) 



where C is constant. The second term in the right-hand side of eq. (22) does not contribute to the 
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amplification, though it enhances and suppresses the amplitude temporarily. For D = 0, 1 have 

(23) 



z + C 



f dzcr(z) = -s (z - zo) + In 



zo + C" 

where C" is constant. The second term in the right-hand side of eq. (23) does not contribute to the 
amplification substantially. 

Next, I consider the case that the quantity r(z) is time dependent. For such the case, the region 
[zo,z] is divided into small regions with the width of time interval Az. Moreover, the region with the 
width Az is divided into quite small N regions numbered '/ in which the quantity r is constant. I 
define the quantity AWj by rjAz/N, where rj is the value of r in the region '/. This quantity AWj is a 
wiener process and the distribution function of AWj is given by 

1 ( (AWj) 2 ) 



P(AW;) = , exp 
3 V2tt(Az)/7V 



vwsY <24) 



N 

Then the quantity AW = ^ AWj obeys the distribution function P(AW) which is given by 

7=1 



1 / (AW) 2 \ 

P(AW) = exp - ^— -L . (25) 



^JlrfM \ 2(Az) 

Therefore, the values of AW in the regions with the width Az are distributed with the probability 
P(AW). Then the statistical average of a variable O is taken as follows: 

y~«00 

(O) - d(AW)P(AW) O. (26) 



From eqs. (16), (21), (22) and (23), the exponent of tp in unit time of z (I denote §) should be estimated 
by 

0= f d(AW) P(AW)@(D) (27) 

where &(x) is the step function which is 1 for x > and for x < , and D = As 1 - 4or(AW)/Az. This 
integration can be performed and I obtain the following expression of Q\ 

a W I (Az)/\ / (Az) l ' 2 s 2 \ 

<* = vn^m *>* j [ —J ' ^ 

where D v is the parabolic cylinder function. 27 ' 28) I note that the expression Q depends on Az. Never- 
theless, I can read the behavior of the exponent as a function of a from this expression. The quantity 
Q as a function of a has an extremum which is determined by dQ/da. I use the following formulas of 
the parabolic cylinder function in order to simplify the expression of dQjda: 

^D p (x) - ^xD p (x) + D p+l (x) - 0, (29a) 
D p+l (x) - xD p (x) + pD p ^(x) = 0. (29b) 
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As a result, I find that the sign of dQ/da depends on only the function Di/2<x). I obtain the subsequent 
condition that § is extremum at a certain a: 



It is known that D v (x) for positive v has [v + 1] zeros, 29 ) where [v + 1] is the maximum integer which 
is not greater than (v + 1). Then the equation, Di/2(x) = 0, has one solution, and I write the solution 
as x so \. This solution is negative, and then a is positive at the extremum of Q. Therefore Q has one 
extremum surely at a positive a. The value at the extremum of the exponent § is given by 



It is shown that Q mia is smaller than s, and then Q has one minimum. This indicates that the exponent 
Q is suppressed by white noise when the value of a is adequate. I must note that the expression Q mm 
is independent of Az. Equation (28) for quite small s should be invalid, because the approximation of 
CD given in eq. (16) does not work well. 

3. Numerical calculation of the exponent by a symplectic method 

In this section, I attempt to solve eqs. (5a) and (5b) numerically. Our purpose is to investigate the 
amplification of the amplitude of when white noise acts multiplicatively. Therefore, the amplitude 
must be calculated precisely, at least, when no periodic coefficient and no white noise term exist. It 
is well-known that the structure of such a system is symplectic. In the similar way, it was proved that 
the system has the symplectic structure even when noise exists if some conditions are satisfied. 30 ) 
Taking this property into account, I use the symplectic method developed in ref . 3 1 in order to solve 
the stochastic differential equations with multiplicative white noise. The first-order method given in 
ref. 31 is applied to the equations in this paper. The initial conditions are 0(0) = 1 and 0(0) = in 
these calculations. The equations are solved numerically from z = to z - 500. The time step in z is 
set to 0.05. 

In the case of a + 0, one trajectory of <p(z) can be calculated when the sequence of noise is given. 
I calculate many trajectories and take their average to obtain the mean value of the trajectories of the 
variable 0- ;) (z), where the subscript i indicates the batch and the superscript (J) indicates the trajectory 
in a certain batch i. In the present calculation, one batch contains 500 trajectories and 20 batches are 
taken. I calculate the mean value A4 ; (z) of the trajectories in the batch i. The mean value over 20 
batches, 0(z), is given by 



I note that it is possible to perform interval estimation by using and Mi. In the case of a - 0, there 




(30) 




(31) 




(32) 
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Fig. 1. Exponents on the y-/3 plane for various values of a. The exponents are calculated by solving the 
stochastic differential equations numerically by the symplectic method. The parameters are (a)a = 0.0, 
(b)ff = 0.5, (c)a = 1.0, (d)a = 1.5, (e)« = 2.0, (f)a = 2.5 respectively. 

is no need to calculate many trajectories. Therefore one trajectory is calculated numerically. 

The exponent is estimated from the average cf>(z) in the range of 200 < z < 500 in order to decrease 
the effects of the initial conditions. This estimation is performed as follows. 1) the set (zk, In (f>(zk)) is 
determined, where z.k is the time at which (f>(zk) is a local maximum and positive. 2) the set is fit with 
a linear function. The coefficient of the time z is adopted as the exponent. Here, I note the reason why 
the values, In 4>(zk), are fit. One way to estimate the parameters is to fit the average <p(zk) directly. In 
such the method, it is implicitly assumed that the dispersion of the distribution of the data at time z, and 
that at time z' (+ z) are the (approximately) same. However, the dispersion is wider with time z in the 
present case, because the treated process is a wiener process. The effects of non-equivalent dispersions 
are decreased by taking the logarithm of the data. Therefore the transformed data, In 4>(zk), are fit with 
the linear function. I may notice that the exponents extracted by the above procedure are different 
generally from the Lyapunov exponents which are estimated by the mean value of the logarithm of 
(f^P . The quantity, ln^(^), is calculated, because I focus on the enhancement of the variable cf> in this 
paper. 
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a a 

Fig. 2. Exponents on the resonance bands. The cross represents the data estimated from numerical results 
obtained by solving the eqs. (5a) and (5b). (a) The values of the parameters, /? and y, are both 2. (b) The 
values of the parameters, /? and y, are 2 and 0.9 respectively. 

Figure 1(a) is the map of the exponents of Mathieu equation, eq. (7), on the y-/3 plane. The step 
sizes in y and /? in the numerical calculations are taken to be 0.02 to draw this figure. I denote these 
step sizes as Ay and A/3 respectively. The color of a square is determined from the arithmetic mean of 
the exponents at four corners which are located at (y,/3), (y + Ay,/?), (y,/? + A/3) and (y + Ay,/? + A/3). 
The resonance band around y = 2 corresponds to the first resonance band of eq. (6). In the same way, 
the resonance band around y = 1 corresponds to the second resonance band. The nth resonance band 
of eq. (6) corresponds to the band around y - 2/n, where n is positive integer. 

Next, I show the map of the exponents for various values of a on the y-/3 plane. Figure 1(b) is the 
map at a - 0.5, (c) is at a - 1.0, (d) is at a = 1.5, (e) is at a - 2.0, and (f) is at a - 2.5. The step sizes 
in /? and y are 0.05 in the numerical calculations for Figs. 1(b), 1(c), 1(d), 1(e) and 1(f). The color of a 
square is determined in the same manner as in Fig. 1(a). The resonance band of the Mathieu equation 
is shown in Fig. 1(a). As shown in Figs. 1(b), 1(c), 1(d), 1(e) and 1(f), the band structure is destroyed 
by white noise, and the values of the exponents become large with a for many sets of (y,/3). However 
it seems from these figures that the exponent on the resonance band is not a monotonically increasing 
function of a. Moreover, the /? dependence of the exponent in Fig. 1(f) is weak as compared with those 
in other figures: Figs. 1(a), 1(b) and 1(c). This implies that the values of the exponents of the equation 
with the periodic coefficient are close to those without the periodic coefficient in Fig. 1(f). (The values 
of the exponents at /? - correspond to the values in the case of no periodic coefficient.) It is evident 
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that the effects of the periodic coefficient become weak relatively. 

Furthermore, I investigate the a dependence of the exponent on the first and the second resonance 
bands. I draw the a dependence of the exponent with fixed parameters, y and ft. I show the exponents 
for the set (y = 2, ft = 2) on the first resonance band, and the set (y = 0.9, ft - 2) on the second 
resonance band. Figure 2 shows the a dependences of the exponents. The cross represents the data 
obtained by solving eqs. (5a) and (5b) numerically. The suppression by white noise is obviously seen 
and there is only one local minimum in each figure. The exponent decreases with a and reaches the 
minimum. It continues to increase with a after that. This behavior is interpreted as follows. The growth 
of the amplitude depends on the mechanism of parametric amplification for small a. This mechanism 
is destroyed by noise with the increase of a. Then the exponent decreases with a. Contrarily, the 
amplitude is amplified by noise for large a, as shown in many researches. In summary, the exponent 
decreases with a, reaches the minimum, and increases after that. The exponents for other parameter 
sets, (y,yS), on the resonance bands behave similarly. 

Finally, the minimum value of the exponent as a function of a is estimated for various values of ft 
and y. I denote the minimum value of the exponent estimated numerically as s m [ n . Obviously s m [ n is a 
function of /3 and y. In this calculation, the range of a is set to [0, 2] and the step size in a is set to 0.01. 
The range of y is set to [0.7, 2.7] and the step size in y is set to 0.5. The exponents for various values 
of a with fixed /3 and y are estimated and s m i n is set to the minimum value of these exponents. I plot 
the quantity s m [ n /s, because the exponent at a = 0, s, is also a function of /? and y. I show the values 
Srma/s for s > 0.3 in order to compare them with the value Q m m/s. The value @ m mls is approximately 
0.893 from eq. (31). 

Figure 3 shows the values Smin/s for s ^ 0-3 and the exponents s. The parameter /? is set to 2.0 in 
Fig. 3(a) and 1.5 in Fig. 3(b). Cross represents data points of s m i D /s and broken line indicates @ m mls. 
Asterisk represents data points of s. As seen in Fig. 2, noise influences the values in the numerical 
calculations. Therefore it is likely that the ratio s m i n /s fluctuates and that the values s m i n /s around the 
maximum of s are below the value Q m i n /s. Then I calculate also the simple moving average of the 
exponents as a function of a with fixed ft and y, and attempt to find the minimum value of them. I take 
the average of n adjoining exponents and denote this average as s™ A ". For an example, the minimum 
of the averages of three adjoining exponents is represented as s™ A3 - Figure 4 displays the s™ A3 /s 
for s > 0.3 and the exponents s. The parameter /? is set to 2.0 in Fig. 4(a) and 1.5 in Fig. 4(b). The 
symbols in Figs. 4(a) and 4(b) are the same as in Figs. 3(a) and 3(b). It is found from Figs. 3 and 4 that 
0mmls is close to the values estimated by numerical calculations around the peaks of s in the resonant 
regions. 
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Fig. 3. The values s m i n /s for various values of y. The parameter f3 is set to (a) ft = 2.0 and (b) /3 = 1 .5. Cross 
represents the data points of s min /s and asterisk represents the data points of s. Broken line is the value 
@ m m/s that is approximately 0.893. 




Fig. 4. The values i SMA3 / s for various values of y. The parameter /? is set to (a) [3 = 2.0 and (b)/3 = 1 .5. Cross 
represents the data points of s m [ n /s and asterisk represents the data points of s. Broken line is the value 
Gmmh that is approximately 0.893. 



4. Discussion and Conclusion 

I studied the amplification in a Mathieu-like equation in the presence of multiplicative white noise. 
An approximate expression of the exponent was derived by introducing the width of time interval on 
parametrically amplified regions where parametric amplification occurs when no noise exists. The ex- 
ponents were calculated by solving the stochastic differential equations numerically by the symplectic 
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numerical method. The intensity of noise and the strength of the coupling between the noise and the 
variable are reflected to the value of the parameter a. The behavior of the exponents as a function of a 
is extracted. 

With regard to the effects of multiplicative white noise on the amplification, it was shown that the 
band structure of the Mathieu equation is destroyed when white noise exists. The resonance structure 
survives for small values of a, and this structure is lost for large values of a. In the previous paper, 32) 
I investigated the amplification in a stochastic differential equation without a periodic coefficient, and 
found that the exponent is a monotone increasing function of a. In contrast, it is found in the present 
study that the exponent as a function of a has one minimum on the parametrically amplified region of 
a = 0. This indicates the suppression of the amplification by white noise which occurs when the value 
of a is adequate. Equation (28) can roughly explain the behavior of the exponent as a function of a: 
The exponent decreases with a, reaches the minimum and increases after that. 

One expects that the exponent as a function of the intensity of white noise has one minimum intu- 
itively. However the exponent may have some minima caused by white noise. Theoretical expression 
given by eq. (30) with the cylinder function indicates that only one minimum exists. This fact is nu- 
merically supported too. It is shown in the previous sections that the exponent as a function of a has 
one minimum theoretically and numerically. 

The minimum value of the exponent as a function of a is estimated from the numerical calcu- 
lations. I calculated the ratio s m [ n /s: the minimum value divided by the exponent at a = 0, s. This 
ratio obtained numerically is in rough agreement with that obtained theoretically around the peaks of 
s on the resonant regions. The minimum value of the exponent is approximately proportional to the 
exponent s. The relative variation is of the order of 90%, as shown in the figures and eq. (31). It seems 
that the variation is small. Nevertheless, the amplitude of the observable is affected strongly after long 
time passes, because this is the variation of the exponent. 

The decrease of Lyapunov exponent by noise was found in the system of an inverted Duffing 
oscillator with noise. 33) The mechanism of the amplification in the present case is different from that 
in the case of the inverted Duffing oscillator, when noise is absent. However, the mechanism of the 
suppression in the present case is surely the same as that in the case of the inverted Duffing oscillator. 
In both cases, the amplification is suppressed by noise which intensity is adequate. The exponent 
decreases with the intensity of noise, and reaches the minimum. After that, the exponent increases 
with the intensity of noise. The decrease of the exponent by white noise implies the possibility of the 
large decrease by colored noise, though the suppression by white noise is weak in the present case. 
Then the system that parametric resonance occurs may be stabilized by colored noise, as found in the 
system of the inverted Duffing oscillator. 
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The expression of the exponent obtained theoretically includes the artificial parameter Az. Then 
the value of a at the minimum of Q depends on Az, while the minimum value of Q is independent of 
Az. The value of Q for large a are out of accord with that obtained theoretically, because the expression 
of Q is calculated in the case that the exponent s is important. These differences may come from the 
fact that the effects of white noise on parametric amplification are not reflected sufficiently. I would 
like to solve these problems in the future study. 

Appendix: The effects of white noise on the exponent of the stochastic differential equation 
with a periodic coefficient 

In this paper, I treated the stochastic differential equation with a periodic coefficient which varies 
sinusoidally. In the same way, I can discuss a stochastic differential equation with coefficients varying 
periodically. I treat the following equation: 

4>(z) + [h(z) + ar(z)] <p(z) = 0, (A- 1) 

where a is constant, r(z) is white noise which properties are given in eq. (4) and the dot represents the 
derivative with respect to z. The function h(z) is a periodic function of period n. 
I treat the subsequent equation to analyze the solution of eq. (A-l): 

O(z) + h(z)®(z) = 0. (A-2) 

From Floquet's theorem, this differential equation has a solution as follows: 

(D( Z ) = e ^ z F(z), F(z) = F(z + n). (A-3) 

Substituting 0(z) = (&(z)if/(z) into equation (A-l), I obtain 

if/ + 2 ($/o) ifj + ar{z)ip = 0. (A-4) 

By using eq. (A-3), eq. (A-4) is given by 

if/ + 2(ifi + F/F) ip + ar(z)if/ = 0. (A-5) 

As shown in §2.2, 1 obtain an approximate differential equation for if/ when the average of the coeffi- 
cient of i// in the period n is taken: 

if/ + 2/yUtA + ar{z)if/ = 0. (A- 6) 

I treat the case that [i is pure imaginary, because I study the case that the amplitude of <E> grows. It 
is conjectured that the effects of the constant coefficient of ip, 2ifi, on the growth of if/ is stronger than 
those of the periodic coefficient of ip, 2F/F. 



13/15 



J. Phys. Soc. Jpn. Full Paper 

Setting \i as -is where s is positive constant, eq. (A-6) is rewritten as 

fi + 2s<j/ + ar(z)4f = 0. (A-7) 

This equation is just eq. (18). Then the final expression of the exponent of <f>(z) is equal to eq. (28). 
Therefore the exponent as a function of a has one minimum on the region where O grows exponen- 
tially. This result suggests that the amplification of the variable <p(z) that obeys the stochastic differen- 
tial equation with the periodic coefficient, eq. (A-l), is suppressed by white noise when the variable 
cD(z) grows and the parameter a is adequate, as shown in §2.2 and 3. 
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